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ABSTRACT 

We show that fuzzy spheres are solutions of Lorentzian IKKT matrix models. The solutions serve 
as toy models of closed noncommutative cosmologies where big bang/crunch singularities appear only 
after taking the commutative limit. The commutative limit of these solutions corresponds to a sphere 
embedded in Minkowski space. This ‘sphere’ has several novel features. The induced metric does 
not agree with the standard metric on the sphere, and moreover, it does not have a fixed signature. 
The curvature computed frorn the induced metric is not constant, has singularities at fixed latitudes 
(not corresponding to the poles) and is negative. Perturbations are rnade about the solutions, and are 
shown to yield a scalar field theory on the sphere in the commutative limit. The scalar field can becorne 
tachyonic for a range of the parameters of the theory. 
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1 Introduction 


It is well known that fuzzy spheres and fuzzy coset spaces[T]-|8] are solutions to matrix models. More 
specifically, they are solutions to the bosonic sector of Euclidean Ishibashi, Kawai, Kitazawa, Tsuchiya 
(IKKT) matrix models.0 An application of these solutions to particle physics has been to rnake extra 
dimensions noncommutative. m Here we show that fuzzy spheres can also be solutions to IKKT matrix 
rnodels with a Minkowski background metric tensor. This means that in addition to making extra 
dimensions noncommutative, fuzzy spheres and more generally fuzzy coset spaces can be used to make 
space-time noncommutative. Moreover, they can be toy models for noncommutative cosmological space- 
times. 

Various aspects of Lorentzian IKKT matrix rnodels have been discussed in the literature, including 
classical solutions and their implications for cosmology. EMiamiiaMisi The solutions were gener- 
ally written in terms of infinite-dimensional matrices, and they may or may not be associated with 
finite dimensional Lie algebras. The advantage of now having fuzzy spheres and fuzzy coset spaces as 
solutions is that they can be expressed in terms of N x N matrices, where N is finite, and they lead to 
closed space-time cosmologies upon taking the N —» oo lirnit, corresponding to the commutative lirnit. 
Big bang/crunch singularities should then appear in this limit, while the finite dimensional matrix 
description is singularity free. 

In section two we write down a fuzzy sphere solution to the Lorentzian IKKT model. The model is 
written down specifically in three space-time dimensions and cubic and quadratic terms are included in 
the action. We show that the solution yields a closed (two-dimensional) universe in the commutative 
limit. While the commutative limit of the solution is topologically a two-sphere, there are a number of 
novel features, arising frorn the fact that it is embedded in a three-dimensional Minkowski space. The 
induced metric does not agree with the standard metric on the sphere, and moreover, it does not have a 
fixed signature. The curvature computed from the induced metric is not constant and it is negative. It 
is singular at two fixed latitudes (which are not located at the poles) and time-like geodesics originate 
and terminate at these latitudes. Thus in this toy rnodel, the big bang/crunch singularities occur at 
nonzero spatial size. 

We examine perturbations around the fuzzy sphere solution in section three. In tlie commutative 
limit, the perturbations are described by a scalar field coupled to a gauge field. The latter can be 
eliminated yielding a scalar field which can propagate on the Lorentzian region of the two-dimensional 
surface. Depending on the choice of parameters the scalar field can be massive, massless or tachyonic. 

Concluding remarks are given in section four. 


2 Fuzzy sphere solution to the Lorentzian IKKT matrix model 

The setting here is the bosonic sector of the Lorentzian IKKT matrix model in three space-time dimen- 
sions. The dynamical degrees of freedom for the matrix model are contained in three infinite-dimensional 
Hermitean matrices /ii = 0,1, 2, with /r = 0 indicating a time-like direction. In addition to tlre stan- 
dard Yang-Mills term, we include a cubic term and a quadratic term in the action (which are both 
necessary for obtaining fuzzy sphere solutions) 

S(X) = ±Tr(-±[Xr,X v ][X*,X v ] + liae lu , x X' t X''X x + £x it X' i ') , (1) 
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where g, a and /3 are real coefficients. Our conventions are £012 = 1, and we raise and lower indices 
/r, v ,... with the flat metric [r/^] =diag(—1,1,1). The resulting equations of motion are 

[[* M , X A,X V } + iae livX [X v , X x ] = -pX^ (2) 

The dynamics is invariant under three-dimensional Lorentz transformations, X M —>• L^ v X v , where L 
is a 3 x 3 Lorentz matrix, and unitary ‘gauge’ transformations, X M —► UX^U', where U is an infinite 
dimensional unitary rnatrix. The equations of motion also have discrete symmetries, namely proper 
reflections. An example is 

(X°,X\X 2 ) -+ (—X°, X 1 , —X 2 ) (3) 

Translation invariance in the three-dimensional Minkowski space is broken when /3^0. 

When /3^0, there exist finite dimensional matrix solutions to the equations of motion © , which 
are associated with tlie su(2) algebra in an N —dimensional representation. Say the latter is spanned 
by N x N hermitean matrices Ji, i = 1,2,3, satisfying [</,;, Jj\ = iaeijk Tfcj^] Let us set 

X° = ^J 3 A' 1 = — Ji X 2 = ^J 2 , (4) 

a a a 

where w, are real. Upon substituting this expression into the equations of motion one gets 

(w\ + w 2 + /3)wz + 2aw\w 2 = 0 

(w 2 — w 2 + (3)wi — 2aw 2 W3 = 0 

(w\ — w\ + (3)w 2 — 2awiw 3 = 0 , (5) 

which has nontrivial solutions. Lorentz symmetry is in general broken by the solutions, unlike the case 
with de Sitter and anti-de Sitter solutions.[ r 14j.[?] The su(2) Casimir operator for any of the solutions 
can be written as -=z(X 0 ) 2 + ^(X 1 ) 2 + Xj(X 2 ) 2 , which has the value \(N 2 — 1) in the N —dimensional 
representation, thereby defining a fuzzy sphere, or actually fuzzy cllipsoid since rotational invariance in 
the (X°, X 1 , X 2 ) space does not in general hold. 

We note that the solution is invariant under the three-dimensional rotation group (rather than the 
Lorentz group) in the special case where w\ = w 2 = w 2 . Let us rnore generally restrict to the case of 
rotational invariance in the (X^X 2 ) plane, which means w\ = w\. Two simple solutions exist in this 
case: 


0 

II 

to 

X 1 - ^ Ji 

x 2 - ^ .h. 

(6) 


a 

a 


and 



0 

II 

to 

x 1 = J, 

a 

x 2 =' /=3 A, 

a 

(7) 


Nontrivial solutions require the presence of both the cubic and quadratic terms in 0, a/0 and (3 < 0. 
@ and (7} are equivalent due to the discrete symmetry ©• For tlie sake of definitess we choose to 
work with the former ([6]). The su( 2) Casimir operator for this solution can be written as 

- ^(X 0 ) 2 + (X 1 ) 2 + (X 2 ) 2 , ( 8 ) 

§The Levi-Civita symbol here is associated with Euclidian space, unlike the ones appearing in Q and ([2]) which are 
associated with Minkowski space. 
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having the value —j(N 2 — 1) in the N —dimensional representation. The ‘time’ matrix X° tlien has 
discrete eigenvalues 2am, where m = ~^ +1 , ~ j ), + 3 ,^=1. For any m defining a time-slice we can 
also define a spatial size. Call A the ‘space’ matrix, where A 2 = X + X_ and X± = X 1 ± iX 2 . We can 
identify it with — =i(J 2 — Jf — J 3 ) for the solution (jçj). A 2 then commutes with A'° and has eigenvalues 
— (3^ N — m 2 — mj. Thus time and the spatial size are discrete. Examples of spectra for (X°,A 2 ) 
for some N —dimensional representations are 


N = 2 

(-«,-£), (a, 0 ) 


N = 3 

(—2a, —2/3), (a,-/3), (2a,0) 


N = 4 

(—3a, —3/3), (—a, —4/3), (a, -3/3), (3a,0) 


N = 5 

(—4a, —4/3), (—2a, —6/3), (a,-6/3), (2a,-4£), (4a,0) 

(9) 


Say a > 0. Then for large N, the spatial size operator A has eigenvalue \J—j3N for the lowest time 
eigenvalue ~ — aN, i.e., the initial state. It then increases to the maximum value yJ—/3 N /2 as the 
time goes to zero, ancl then decreases to zero upon approaching at highest time eigenvalue ~ aN, 
i.e., the final state. This solution can thus be regarded as a discrete analogue of a closed cosmological 
space-time. The eigenvalues of X° versus those of A are plotted for N = 100 in figure 1. 



Figure 1: Fuzzy closed universe solution. Plot of the eigenvalues x n of the tirne matrix X° versus tlre 
eigenvalues a n of the space matrix A for N = 100, a = 1 and /3 = —1. 


Just as with the fuzzy sphere in a Euclidean background, the commutative limit of the matrix 
solution here is obtained by taking N —>■ 00 . Here we also need a,/3 —► 0, with aN and \J—(3 N finite 
in the limit. The commutative limit of the solution is tlren characterized by two real parameters, which 
we denote by oq and r 2 , 


—> r 


( 10 ) 


One standardly defines the commutative limit in analogous fashion to the classical limit of a quantum 
theory, where we can take a to play the analogous role to Ti. This means replacing the matrices X^ by 
space-time commuting coordinates which we denote by a±, where x° and x l , i = 1,2 denote the time 
and space coordinates, respectively. The constraint on the su( 2) Casimir operator ([ 8 ]) means that in 
the commutative limit the solution satisfies. 


a 2 0 (xy + (xr + (xy = 


( 11 ) 
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While real clq means that the solution is topologically a two sphere, there are a nurnber of novel features, 
which we show below, due to the fact tliat tliis ‘sphere’ is embedded in Minkowksi space-time. 

The commutative limit also requires replacing the commutator of functions of , evaluated for the 
solution ([6]), by ia times Poisson bracket of the same functions of the coordinates a+. Thus the Poisson 
brackets of the embedding coordinates are 

{x°,x x } = -2x 2 {x 2 , a; 0 } = — 2X 1 {x 1 ,x 2 } =—2agX° (12) 


We can express a+ in terms of angular momenta ji, 
algebra {ji,jj} = djkjk, using 



* = 1,2,3, which satisfies the su{ 2) Poisson bracket 

/ h \ 

I a oji I , (13) 

V — «oj '2 / 


and from (111 , j 2 + j 2 + j\ = (t^) 2 . l" 01 simplicity, we set r = 2ao so that j t spans a sphere of unit 


radius. We can introduce standard spherical coordinates {9, (j>), 0 < 9 < n, 0 < <j> < 2tt, and write 


j i = sin 9 cos (f) j ‘2 = sin 9 sin <j> j% = cos 9 


(14) 


The su{ 2) Poisson bracket algebra for j, is recovered upon defining the Poisson brackets on the sphere 
to be 

{F, G}{9, <t>) = csc 9 (deFd+G - d^FdgG) , (15) 

for any two functions F and G on the sphere. 

The induced metric g ab = daX^dbX^, a,b ,... = 9, <j>, does not agree with the standard metric on tlre 
sphere, and moreover, it does not have a fixed signature. The curvature computed from the induced 
rnetric is not constant and is negative! The invariant interval constructed from the induced metric is 

— dr 2 = 4^ag cos 2 9 — sin 2 9^j d9 2 + 4ag sin 2 9 d(f> 2 (16) 

ggg vanishes at two latitudes 9 = 6± on the sphere defined by tan0± = |ag|. Say that 9 = 9± 
is contained in the northern hemisphere, 0 < 9 + < j, while 9 = 9- is contained in the southern 
hemisphere, ± < 9 < ir. The signature on sphere is Euclidean for 0 < 9 < 9 + and 9_ < 9 < £, while 
it is Lorentzian for 9 + < 9 < 9_. We can regard 9 as a time-like variable for the latter, with 2ao sin6* 
being the spatial radius at any time-slice. 9 = 9± correspond to singularities in the curvature, and are 
not coordinate singularities. The Ricci scalar computed frorn the induced metric is 


R = - 


2(oq cos 9 2 — sin 2 9) 2 


(17) 


and thus it is singular at the latitudes 9 = 9±. ( |17| shows that the curvature in the nonsingular 
regions is everywhere negative. The singularities of the Ricci tensor are analogous to big bang/crunch 

2a 2 

singularities, with the distinction that they occur at a nonzero spatial radius 2ag sin0± = / on the 

V a o + 1 

two-dimensional space-time. Time-like longitudinal geodesics originate and terminate at the singular 
latitudes 9 = 9±. This is because their tangent vectors (^|, ±p) = ( , . 2 1 2 2 , 0) are well 

defined in the Lorentzian region, 9± < 9 < 9_, while they are imaginary in the Euclidean regions, 
0 < 9 < 9± and 9- < 9 < |. The total elapsed proper time along these geodesics is finite and given by 

the elliptic integral 2 / t a n -i” 0 a ° d9^J sin 2 9 — cos 9 2 . 
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3 Emergent field dynamics 


Here we perturb around the matrix solution ©■ Similar to [S], we can define noncommutative field 


fiu 


Here we take 



II 

O 

V 0 ,* 1 ] 

a 

+ 2iX 2 


O 

to 

II 

~[X°,X 2 ] 

- 2iX x 


F 12 = 

V[A\.Y 2 ] 

a 

- 'Jx° , 

2a 

(18) 


which transform covariantly under unitary gauge transformations, F^ v —► UF^U', and vanish when 
evaluated on the fuzzy sphere solutions §• The matrix action ([!]) can then be re-expressed in terms 
of the noncommutative field strengths 


g 2 S(X) = Tr {—F^F^ - -ia 2 (F m X 2 + F^A 1 ) + ia 2 [\ - ^=)f 12 X° 
^43 V3 2 or/ 

+ (f-^)((x 1 ^ + (x V ) + ^-5) ( xT} (i 9 ) 

Now perturb around the matrix solution § using 

A'° = 2 (j 3 + -^=A°) X 1 = Ji _). aA 1 X 2 = j 2 _ a A 2 , (20) 

V \J—p ) a a 


where the perturbations are functions on the fuzzy sphere, = A^(Ji, J 2 , J 3 ). If we write infinitesimal 
unitary gauge transformations using U = 11-t=A, where A is a hermitean matrix with infinitesimal 




elements, then the infinitesimal variations of read 

8A» = -t(i[A,J'*] + -^=[A,A'‘]) , 


where we identify (J°, J 1 , J 2 ) with (J 3 , Ji, J 2 ). Substituting (20) into (19) gives 

S(X) = “jTr| - l ^ v F» v -^ia(F m A 2 +F 20 A 1 )+‘^Ç--^')F l2 A° 


8 ia 
”3 

'P 2a 2 


U+ ([J 1 .+ 2 ] - [4,+‘])-4“ - 2 io(j - A) [j4 i,aV 3 


(f - ^) ((+‘) 2 + (+ 2 ) 2 ) - 2a 2 (jA _ 5)(X»)»] + S(X| 


solution/ 


( 21 ) 


( 22 ) 


As stated previously, tlre commutative limit is obtained by taking N —> 00 , along with a, /3 —> 0 and 
both aN and \J—fiN are finite in the limit. Upon using ( 10 ) and (13), the commutative limit of the 
field strengths ([l 8 | is 


F m -+ ^({^A 1 }-^,^ 0 }-^ 2 ) 
F 02 -+ -2ia({j 5 ,A 2 }-{j 2 ,A°} + A^ 
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(23) 


F 


12 


-2iaa 0 ({ji,A 2 } - {j 2 , -4 1 } - -4 0 ) , 


where are now functions on the commutative sphere. The trace on functions of the fuzzy sphere 
is replaced by the corresponding integration on the sphere in the commutative lirnit. The relevant 
integration measure d/i(0, <fi) slrould be such that tlie standard trace identities survive in the lirnit, i.e., 
for any three functions G,H and I\ on the sphere we want f (j){G, H}K = f dfi(6,(f))G{H,K}. 
From (|~i~5j) we need to choose the standard integration measure on the sphere dfi(6, <j>) = sin ddddcj) 


(rather than say ,J—gdOd<f), where g is the determinant of the induced metric). Then the action (22) 
reduces to 


S(X) — S(X | so l u ti on ) —> —=- 


2 a 

9c 


J sin6d6d<p ^ - ({j 3 , A 1 } - {ji,A 0 }) - ({j 3 , A 2 } - {j 2 ,^°}) 
+ a o({jl,^4 2 } — |j2, ^ 1 }) + 2 (Oq + 1){j 3 , A^}A 2 


+(A 0 ) 2 -, 


j((A 4 ) 2 + (A 2 ) 2 )}, 


(24) 


where g c is the commutative limit of the constant g. Following [Sj we write the perturbations A^ in 
terms of gauge potentials (Ag,A$) and a scalar field j) on the sphere using 

A° = A,p + j 3 ip 


A 1 = —sin<pAg — cot6cos(pA c f, + ji'i/) 
A 2 = cos <pAe — cot 9 sin <pA<f, + jzip 


(25) 


Then from tlre fundamental Poisson bracket (15), gauge variations (SAg,SA,p) = (dgA,d$A) agree with 
the commutative limit of (21), where A is now an infinitesimal function on the commutative sphere. 
Substituting (25) in (24) gives 


S(X) - S(X| solution ) 


f smddddc/) } (oq cot 2 9 — 1)— csc 2 O^d^ip) 2 
9c J \ 

+ (oq sin 2 0 — cos 2 d) (dgip) 2 — (3 — 2(ciq + 1) sin 2 o'jtp 2 

+2 csc 0 ((oq + 1) sin 2 0 — 2a 0 + l) Fg^tp — 2 cos 0(a o + 1 ^Tg^dgip 


(26) 

where Tgp = dgAp — dpAg is the 17(1) gauge field on the surface. We remark that the gauge field and 
scalar field kinetic energies can have opposite signs, a feature that was present in similar two-dimensional 
systems. na.m However, gauge fields are nondynamical in two-dimensions. We can solve for Tgp from 
the field equations, yielding 


FQ(f, — 


cos 0(a o + 1 )dgip — ((dg + 1) sin 2 0 — 2a 0 + l) csc 0 tp 


Oq cot 2 0 — 1 


constant 


( 27 ) 


and substitute back into the action. Upon setting the constant equal to zero, we get 

(dgip) 2 csc 2 0 


S(X) - 5(X| solution ) 


2a 4 ag 

9c 


sin 0d0d<p 


(d^ip) 2 - 4 ml ñ ip 2 


{ (a 0 + 1) sin 2 0 — a 0 
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(28) 


16 a 4 aQ 


J sin6»d6»#|-ia a V'9 a V> - |m^ 2 J , 


where the index a = (0, <p) is raised and lowered using the induced metric given in (16). The effective 
mass squared of the scalar field is 6 —dependent 


“eff = 


(oq — 1 )^(a n + 1) sin 2 0 — 3 d\ 
4a„ (J(ap + 1) sin 2 9 - a n j 


(29) 


As stated before, the signature of the induced metric is Euclidean when sin 2 9 < 
when sin 2 8 > 


2 ° i , and Lorentzian 

a n i 


a n + 1 ' 


. Therefore (28) describes a Euclidean field theory for the former and a Lorentzian 


field theory for the latter. There are three different possibilities for the Lorentzian field theory: 


a) Tlie action describes a tachyon when a^ > 1. Tliis is since the factor (ajj + 1) sin 9 — 3oq in (29) is 
negative in this case. 

b) The scalar field is massless when a^ = 1. 

c) The effective mass-squared for the scalar field is positive when 


a n < 1 


ancl 


a n + 1 


< sm 8 < 


3Qq 

a n + 1 


(30) 


It follows that the action (28) describes a massive scalar field throughout the entire Lorentzian region 
when | < a n < 1. On the other hand, when a n < / the scalar fielcl becomes tachyonic in the region 
where sin 2 9 > . 


4 Concluding Remarks 


We found fuzzy sphere solutions to the Lorentzian IKKT model which provide toy models of a non- 
commutative two-dimensional closed universe, where time and spatial size have discrete values. Sin- 
gularities in the Ricci tensor appear in the large N (i.e., commutative) limit. They are analogous to 
big bang/crunch singularities, with the novel feature that they occur at nonzero spatial size. Perturba- 
tions around the fuzzy sphere solution are described by a scalar field in the commutative lirnit which 
can propagate in the Lorentzian region of the manifold. The scalar field can be massive, massless or 
tachyonic, the choice depending on the parameter a n (and also on the range of 9 when a n < 2 ). For 
^ < a n < 1 the scalar field is always massive, ensuring the stability of the commutative field theory in 
this case. Corrections to the commutative limit are obtained by expressing the rnatrix product in the 
action (22) in terms the star product on the sphere|3]-[7| and keeping the next order terms in the 1 /N 
expansion. 


For a more realistic model of a noncommutative cosmological space-time, one can look for fuzzy 
coset space solutions to the IKKT matrix rnodel associated with dimension d > 2.(7 One possible 
example worth consideration is the fuzzy analogue of the four-dimensional coset SU(3)/U(2). For coset 
spaces with d > 4 one may be able to rnake both four-dimensional space-time and extra dimensions 
noncommutative. Just as with the example of the fuzzy spliere, the commutative limit may lead to 
a manifold divided up into regions with different signatures of the metric. Perturbations about such 
solutions are expected to be described by a coupled gauge-scalar theory in the commutative limit A 
common feature of the emergent field theories in previous examples [T5! is that scalar field and gauge field 
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kinetic energies can appear with opposite sign, which also can be seen in (26 1 . This sign discrepancy 


was harmless for d = 2, since the gauge field could be eliminated. On the other hand, it is of concern for 
d > 2 , so it would be interesting to see if this discrepancy can be cured upon taking the commutative 
limit of higher dimensional fuzzy coset space solutions. 
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